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Abstract 

Is it actually possible to interpret gravitation as space's property in a pure classi- 
cal way. Then, we note that extended self-gravitating system equilibrium depends 
directly on the number of dimension of the space in which it evolves. Given those 
precisions, we review the principal thermodynamical knowledge in the context of 
classical gravity with arbitrary dimension of space. Stability analyses for bounded 
3D systems, namely the Antonov instability paradigm, are then rapproched to some 
amazing properties of globular clusters and galaxies. 
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1 Nature of classical gravitation 

1.1 Gravitational field 



In a classical way, gravitation is a force F = — gradf/ acting on a test mass m, 
deriving from a scalar potential energy U which at any time t depends only 
on the position r G M'^. This potential is generated by a mass density field 
p (r, t) such that Poisson equation holds 

U{r,t) = mip{r,t) and A^^p {r,t) = AnGp {r,t) (1) 

Inverting the laplacian A3, one writes the gravitational potential as 

^(r,t) = -G/0^rfr' (2) 
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Using a convolution product, this last relation writes 

ijj ^Att G g3* p 
Where gs is a solution of the equation 

A3 gs^S , 



(3) 



(4) 



which is different from harmonic polynomials. In other words g^ is the laplacian 
Green's function in 

Using information theory, interpretation of equation (3) is clear : Green's 
function is interpretable as impulsional response of the corresponding oper- 
ator, gravitational potential ip (r, t) is then the response from space when it 
is submitted to the presence of a mass density distribution p{r,t). Constant 
factors are also clearly interpretable : the Newton-Cavendish constant G fixes 
units and Att is the value of the surface of unit radius sphere in R^, it is then 
attached to the radial nature of the gravitationnal interaction. 

Let us note that the traditional link usually made between gravity and space 
in general relativity is already contained, using this formulation, in classi- 
cal field gravity. Two main problems remain: the instantaneous character of 
the interaction and the fundamental restriction to isotropic spaces (laplacians 
acts equally in all directions). Two major fulfillments of Einstein theory of 
gravitation were to solve these problems. 

In a classical context, one can directly generalize equation (3) to obtain the 
classical definition of the gravitational potential in a n— dimensionnal space, 
we have 

Ip^G Sn-l gn*p (6) 



where 



Sn-i = — — 7:^ with e 
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r{n/2) 

represents the value of the surface of unit radius sphere in R", and 
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if n = 2 
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if n > 2 
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is Green's function of the lapacian operator in M". 

Representation of the gravitational interaction in M" via relation (6) is not 
original but it is presented here in a rational way. 



1.2 Equilibrium 



Let us describe global dynamical properties of a system of particles of 
same^ mass m represented by their positions and impulsions merged in the 
vector r^^j ^ := (rj(t), Pi{t) = mii). Three tensors are of interest : The kinetic 
/C, potential U and inertial X one, which components are respectively 

Vz, j = 1, ■ ■ ■ , ICij := Uij := mGSn-i^gniji) lij := mViTj . 

3 

In the continuous case where N ^ oo, variable is no longer Fj but its proba- 
bility density / = f(Ti, ■ ■ ■ , r^r, t). Dynamical tensors become usually 



^^3 ■■= 1^ f dr^ := -m I f dr^ 



2m J drj 

r N 

Xij := m TiTj f dT'^ where dF'^ = H 

J i=i 



N 

dVi 



In the general conservative case / obeys Liouville equation, which reduces 
under generic symetry assuptions to CoUisionless Boltzmann Equation for 
gravitating systems (see for example). 

It is not too long so, to prove ^ the fundamental virial theorem 

Theorem 1 If U is homogeneous of degree a, i.e. VA G M, U{Xri, ■ ■ ■ , Arjv) = 
A"W(ri, ■ ■ ■ , tn) then 



It is quite natural to define an equilibrium state by 

d^Ti (X) 



dt^ 







^ this assumption is not essential but it simplifies notably the notations, mainly in 
the continuous case where N ^ oo 

^ from Newton fundamental principle in the discrete case or from Liouville equation 
for the continuous case ... 
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in some mean sense. Hence, for a self-gravitating system in a n— dimensional 
space : 

• li n — 2,U is not homogeneous : Virial theorem does not apply in this form; 

• If n 7^ 2. (/„ then U is an homogeneous function of degree {n — 2) and 
equilibrium is characterized by the relation 



For extended self-gravitating systems like globular clusters or galaxies a 3 
dimensional space in virial theorem seems compatible with observations. 



2 Thermodynamics 

If all particles have the same probability law, are independant and do not 
interact by pair but only globally through their whole mean gravitating field, 
one can reduce the dimensionality of the phase space : The system is sta- 
tistically equivalent to a test particle of mass m at position r G M", with 
impulsion p G M", described at any time t by a distribution function /(r, p, t) 
and evolving in a mean field j/j {r,t) . These two functions are solutions of the 
Vlasov-Poisson system 



2Tr(/C) + (n - 2)Tr(W) = 





2. 1 Definitions 



Several quantities are in general, used in thermodynamics : 



• Phase space variable T 



V = (r, p) G X 



n 



• Space {v) or mass (p) density 




• "Number of particle" N or system's mass M 



N[f] j /(r,i)dr = j v{v,t)dv :=M/ 



m 
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Energy 



E[f] := K[f] + U[f] with 



U[f] :=^/Kr,tMr,t) dr 



Angular Momentum^ 



J[f]- jrp4>nT,t) dr 



• Boltzmann Entropy 



S[f] :^-J f{T,t)ln[f{T,t)] dT 



For precise considerations let us define some set : 
(1) Unbounded Systems 



Gn{N, E) 



f s. t. vr 



K < oo , U < oo 
< oo , < oo 



gn{N,E,j) = g{N,E) n {/ s. t. vr j < 00} 

(2) Bounded Systems Dei?" 

We designe by Supp(/) the support of the distribution function / , i.e. 
the complementary of the largest open set of F e (R" x R") such that 
/ (F) = 0. We then call 

gn{D,N,E) = {fegn{N,E) ; Supp(/) C 
gn{D,N,E,J) = {fegr,iN,E,J) ; Supp(/) C £>} 



2.2 Thermodynamical Equilibrium Problem 



The classical thermodynamical equilibrium problem concerns the existence 
of an entropy extremalizcr. Considering a set ^„ of acceptable distribution 
functions, it can be posed as : 

3? /+ e gn s.t. V/ e gn S[f] < 

^ the quantity indicates the tangential components of the impulsion p, the quan- 
tity r stands for the euchdian norm of the position r of the test particle. 
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For sets considered in the previous section, this problem corresponds to a 
classical Euler-Lagrange one, which solutions are the well-known isothermal 
spheres : 



For / e Gn{E,N) 



= exp I —a — j3(— — h mip' 
2m 



Lagrange multiplier a and (3 correspond respectively to the conservation of 
N andE . 

For / G gn{E,N,J) 

^ =^"^{-"-^(2;^ + ^ + "^^ ^)] 

The additional constraint of J conservation corresponds to the introduction 
of the uj multiplier. 



2.3 2D Thermodynamics : Q2 

We reproduce here the classical results obtained in and The main one 
is twofold : We first produce a bound for entropy, we then study the existence 
of a distribution function which allow reaching that bound. Concerning the 
upper bound for the entropy, we prove that 

Theorem 2 V/ G G2{E,N) 

2E 

S\f]<S+{N,E):= sup S[f]<— + N\n{e7i^) 

Q2{N,E) 

In 2D, entropy of unbounded self-gravitating systems for which E = est, and 
= est is bounded from bellow. 

Theorem 3 In the notations of Theorem 2 

S'^{N,E,J) = S+{N,E) 



Adding the angular momentum's constraint does not change the least upper 
bound on the entropy. The existence of a distribution function /"*" correspond- 
ing to this entropy maximizer, i.e. S [/"*"] = is closely related to the set 
under consideration : 

in the units of original papers 
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In Q2{E, N) : /+ exists and is unique, 



^2 ^g2{_B-Af2)/Ar2 _|_ 

it generates the potential 

V'+ = iVln(e2(^-^')/^'+r2) . 



in the units of papers (|lj) and 
There is no /+ e ^2(^, N, J) for which = S+ . 

In Q2{D, E, N) and ^2(-D, -E, iV, J) : exists and is unique 



2.4 3-D Thermodynamics : Q3 



Introducing a new dimension changes drastically the situation concerning the 
classical equilibrium problem of gravitational thermodynamics. As a matter 
of fact, it is well known for a long time (see for example (5) p. 268) that 

• Entropy has no global maximum on Qz{E, N), Q-i{E, N, J), GsiD, E, N) and 
QsiD^E^N, J). 

• Entropy has no local maximum on Qz{E, N) and QsiE, N, J). 

The existence of local maximum for entropy in Q^{D,E,N) corresponds to 
an extensive literature initiated by the works of V.A. Antonov (see (0) for a 
review) in the early 60's. It represents a beautifull problem of thermodynamics. 



2.4-1 Entropy extremalizer in Q3{D,E,N) 

We reproduce below the main results obtained by T. Padmanabhan (see |3)) 
which clarify all previous works. If we denote by R the radius of the largest 
bowl contained in the spatial part of D. It is proven that any entropy ex- 
tremalizer in QslD, E, N) must be of the form 

N -3/2 

/+ = ( — j Uo e-^^ with muo = p(0) e^^(°) . (7) 

The associated potential verifies Poisson equation 

A3 ^ = AirGmUo e"^'^ (8) 
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with the hmit condition ip (R) = —GM/R. Introducing dimensionless variables 

Lo = ^AnGp (0) (3 Mo = Anp (0) Ll iIj^ = = ^ 
x = r/Lo n = p/p{0) p = M{r)/Mo y = p {^j - (0)) 
Poisson equation becomes 

Milne's functions, v = p/x and u = nx^/p, transform Poisson equation into 



udv 1 — u 
V du M + f — 3 



with 



V = when m = 3 



and 
dv 
du 



-5/3 



{u,v)={3,0) 



Hence, Isothermal extremal spheres lie on a curve in the u — v plane. This 
curve is plotted on Figure 1. 
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Fig. 1. Isothermal sphere in the Milne plane 



2.4-2 Antonov Instability 

In a meaningfuU remark, Patmanabhan notes that dimensionless quantity 

RE If 3\ 
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lies also on the same u — v plane. He then asks the fundamental question : 
Can E^R and M be accomodated by a suitable choice of p(0) and /3 ? As 
one can see on Figure 2, the answer is clearly no. There exists a critical value 
Ac — —0.335, associated with the possibility to put a given isothermal extremal 
sphere in a given box ! 




• If A < Ac isothermal sphere cannot exist, entropy extremum cannot exist; 

• If A > Ac isothermal sphere exist, entropy extremum exist ! 



More exciting is the fact that the nature of extremum depends on k = 
p (0) /p {R) the density concentration ratio of the isothermal sphere. 



• If K > 709 : The extremum of the entropy is an unstable saddle point; 

• If K < 709 : The extremum of the entropy is local maximum. 



This last point is generally associated to the so-called gravothermal catastro- 
phe, we prefer to call it Antonov Instability. As we will see in the next section, 
such an instability is certainly at the origin of some important characteristics 
of extended self gravitating systems. 
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3 Antonov instability in astrophysics 



3.1 Globular clusters in galaxies 



Since early 80's observations have shown that galactic globular clusters split in 
two categories which differ by some properties of their radial density profiles. 
On the one hand a large familly of about 120 clusters with a large constant 
density core which extends to almost the half mass radius (-R50) of the whole 
system. This large core is surrounded by a power law density decreasing halo. 
On the other hand, a small familly of about 20 core collapsed clusters with a 
very high central density which decreases monotonically outward with mainly 
two power law indexes. These two types of globular clusters are very well 
represented by two of their components, namely NGC 6388 for core halo cluster 
and Trz2 for core collapsed cluster, see Figure 1 of 



Such a behaviour of the radial density could be explained in a very simple 
way by Antonov Instability. As a matter of fact, if globular cluster formation 
results from the collapse of a small, hence homogeneous, region of some galaxy, 
the natural result is roughly an isothermal sphere (see (jj)) with generally a 
contrast density k less than the critical value. The evolution of the cluster in 
the galaxy produces a slow evaporation of the cluster (passing through the 
galalactic plane in spiral galaxies for example). Due to the negative specific 
heat of such gravitating systems, this evaporation makes the contrast density 
growing. When k reaches the critical value, Antonov instability triggers and 
transforms the core halo density profile into a collapsed core one. On Figure 3 
we represent radial density profile after collapse of an initially homogenenous 
system and of an initially inhomogeneous one (see (j4)). 



P 



Inhomogeneous 




Fig. 3. Radial density profile obtained from the gravitational collapse of an homo- 
geneous set of mass and an inhomogeneous one 
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3.2 A paradigm for Super Massive Black Hole(SMBH) formation 

From the accumulation of observational data, it becomes necessary to put a 
Super Massive (from 10^ to perhaps 10^ Mq) Black Hole in the dynamical 
center of galaxies. Excepted the fact that such gravitational monsters are 
as older as their host, little is known about their formation process. In the 
context of hierachical galaxy formation scenario, Antonov instability could 
produce a good paradigm. As a matter of fact, if galaxies are the net result of 
successive collapse and merger of gravitational structures, it can be modehzed 
generically by a general collapse of inhomogeneous media. As showed by (4), 
in such case, reaching the center small structures first collapse to form quasi- 
isothermal sphere surrounded by the rest of not yet collapsed large structures. 
Evaporating the smallest, largest's collapse could trigger Antonov instability. 
Progenitor of SMBH could then be formed. This process is up to now to be 
confirmed but seems correspond to all observed properties. 



4 Conclusion 

Classical gravity is an amazing topic. Although, 2D gravitating systems are 
well described by thermodynamics, their equilibrium is not well defined. By 
opposition, provided that n > 3, we possess a powerfull tool to describe gravi- 
tational equilibrium for systems in M", but the corresponding thermodynamics 
is not too efficient. However, in the restricted case of bounded systems in M^, 
the message from gravitational thermodynamics and particulary Antonov in- 
stability, could be fundamental to explain some features of self gravitating 
systems. 



References 

[1] Aly J.-J.,Phys. Rev. E 49, p. 3771, 1994 

[2] Aly J.-J.& Perez J., Phys. Rev. E 60, p. 5185, 1999 

[3] Padmanabhan T., ApJSS, 71, p. 651, 1989 

[4] Roy F.& Perez J., MNRAS, 348, p. 62, 2004 

[5] J. Binney & S. Tremaine, Galactic dynamics, Princeton University Press, 
1987 

[6] Djorgovski, S.; King, I. R.; Vuosalo, C; Oren, A.; Penner, H., Proceedings 
of the 118th. Symposium of the International Astronomical Union, p. 281, 
1986 



11 



